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We propose a new and robust method to test the consistency of the cosmic evolution given by a 
cosmological model. It is realized by comparing the combined quantity rJ MB /Dv N , which is derived 
from the comoving sound horizon ra from cosmic microwave background (CMB) measurements 
and the effective distance Dy derived from low-redshift Type-la supernovae (SNe la) data, with 
direct and independent ra/Dy obtained by baryon acoustic oscillation (BAO) measurements at 
median redshifts. We apply this test method for the ACDM and wCDM models, and investigate the 
consistency of the derived value of rd/Dy from Planck 2015 and the SN la data sets of Union2.1 
and JLA (z < 1.5), and the rd/Dy directly given by BAO data from six-degree-field galaxy survey 
(6dFGS), Sloan Digital Sky Survey Data Release 7 Main Galaxy Survey (SDSS-DR7 MGS), DR11 
of SDSS-III, WiggleZ and Lya forecast surveys from Baryon Oscillation Spectroscopic Data (BOSS) 
DR-11 over 0.1 < z < 2.36. We find that r^ MB /Dv N for both non-flat ACDM and fiat wCDM models 
with Union2.1 and JLA data are well consistent with the BAO and CMB measurements within la 
CL. Future surveys will further tight up the constraints significantly, and provide stronger test on 
the consistency. 


I. INTRODUCTION 

In the last two decades, cosmologists have been making 
great efforts for establishing standard cosmological model 
to describe the contents and evolution of the Universe. 
With precise measurements of the cosmic microwave 
background from the Wilkinson Microwave Anisotropy 
Probe ( WMAP) 0,1 and Planck 11 satellites, the 
cosmological parameters have been measured in a higher 
and higher precision, making it possible to test whether 
the standard A cold dark matter (ACDM) model can de¬ 
scribe the cosmic evolution throughout the history of the 
Universe. In fact, after Planck has published its 2013 re¬ 
sults 1], it was found that the previously measured value 
of Hubble constant Hq through 600 Cepheid variables Q 
is higher than the Planck measured value by 3er confi¬ 
dence level, although later it is shown that by correcting 
the NGC 4258 distance one can obtain a lower value of 
Hq which is compatible with Planck 2013 results §)j. In 
addition, it is shown that the Planck constrained 
parameter plane is in tension with CFHTLenS data 0, 
thermal Sunyaev-Zeldovich effect 0, and statistics of 
cluster number counts [§]. These interesting tensions be¬ 
tween cosmological data sets drive us to consider more, 
and robust test on ACDM model at different periods of 
cosmic evolution. 

A useful and interesting data set of cosmic distance 
estimator is the baryon acoustic oscillation (BAO) data 
from galaxy surveys [jjHS. The measurement of BAO 
scale is normally written as rd/Dy{z) (or Dy(z)/rd). 
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Here ra is the comoving sound horizon at the end of the 
baryon drag epoch, which is completely determined by 
the physics in the early Universe at redshift z > 1100. 
In this regime, the observation of the cosmic microwave 
background radiation can faithfully reflect the physics 
that prevails the Universe. After the photon decoupled, 
the baryons began to fall into the initial gravitational 
potential provided by dark matter, so that the galaxies 
formed in the potential well at the late-time cosmic evo¬ 
lution. Therefore, the gravitational clustering scale of the 
BAO seen in galaxy redshift surveys (~ 105 /i _1 Mpc) is 
fixed and determined by rj and does not depend on late¬ 
time cosmic evolution. On the other hand, the denom¬ 
inator of the BAO measurement, Dy(z), is determined 
by the late-time cosmic expansion, i.e. angular diameter 
distance and Hubble parameter, which can be strongly 
affected by dark energy or modified gravity effect. 

Therefore we would like to perform a test of cosmic 
evolution in a different perspective than Refs. 01]. If 
a cosmology model is the true model of the Universe, 
it should be able to fit observations throughout all red¬ 
shift ranges. Therefore, since Type-la supernova (SN la) 
data is mainly obtained at low redshift (z < 1.5), and 
CMB data at high redshift (z ~ 1100), we will use CMB 
to determine the ra only, and use SN la data to deter¬ 
mine Dy(z). Then we compare the combined quantity 
r^ MB /Z)y N with the direct measurements of rd/Dy(z) 
from BAO surveys at median redshifts out to z ~ 2.4. In 
this way, we calculate the lever arm at two ends (CMB 
and SN la) and compare their prediction in the inter¬ 
mediate redshift regime (BAO). We can then obtain a 
robust test on cosmic evolution by independently consid¬ 
ering different data sets. 

This paper is organized as follows: in Section[Bl we dis- 
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cuss the observational data we use, and show the details 
of the calculation in the comparison for each data set, es¬ 
pecially for the SN la data; in Section ITTT1 we talk about 
the comparison results of r^/Dy and the other quantities 
used in the comparison for the ACDM and u>CDM mod¬ 
els with all three kinds of measurements; the summary 
and conclusion will be presented in the last section. 


II. OBSERVATIONAL DATA 

In this section, we discuss the data we use in the com¬ 
parison, including the SN la data from Union2.1 and 
JLA, the BAO data from galaxy and Lya forest surveys 
and the sound horizon data from the results of Planck 
2015. We also show the details of the derivation of the 
volume weighted effective distance Dy from low-redshift 
SNe la, r^/Dy from mid-redshift BAO and the sound 
horizon 7"d from high-redshift CMB measurements. 


A. SN la data 

There have been a long history of using SN la data 
to constrain the evolution of the Universe. However, in 
recent years, as the constraints from CMB measurement 
have been greatly improved, it is found that some SN 
la samples provide ~ 2er inconsistent results than BAO 
and CMB measurement. For instance, in jjjj, it is found 
that by using the SNLS combined Type-la supernovae 
samples, the joint constraints on wo~w a parameters pre¬ 
fer a phantom dark energy (w < —1) which is incon¬ 
sistent with the Planck+ BAO measurement at 2a CL. 
Same problem is also found in [l3[ and [lij . Interest¬ 
ingly, by re-calibrating the light curve fitting parameters 
with BAO data set, Ref. [l5[ finds that the results of the 
SN constraints shift back to be consistent with Planck 
measurement of the CMB. Thus, there is a possibility 
that the light curve parameters, or the distance-modulus 
relation in SNLS samples have some systematic bias so 
that the results deviate from the true values at 2er CL. 
To avoid or reduce this possible bias, in our approach 
we use the SN la data from Union2.1 catalog 1 [lg, [l7| 
and the Joint Light-curve Analysis (JLA) SN la data 2 
0. These two data sets contain the largest SN la sam¬ 
ples at present, and could provide stringent constraints 
on the cosmological parameters and the quantities we are 
interested in. 

The Union2.1 sample contains 580 SNe la which are 
calibrated by SALT2 light-curve fitter |19|, and its red- 
shift range covers 0.015 < z < 1.414 jl6j]. Following 0, 
the observed distance modulus of Union2.1 data set is 


1 http://supernova.lbl.gov/union 

2 http://supernovae.in2p3.fr/sdss_snls_jla/ReadMe.html 


given by 

^Union = m max + a . Xl _ 0 . c + ■ P - M B , ( 1 ) 

where m^ ax , X\ and c are the three parameters of light 
curve that are fitted by SALT2 [Tij |. The m]§ ax is the 
rest-frame R-band peak magnitude, x\ and c are the 
light-curve shape and color parameters respectively. The 
parameter P denotes the possibility that SNe la belong 
to the host galaxy with mass less than 10 10 TOq, which 
takes account of the correlation between the SN la lu¬ 
minosity and host galaxy mass. The Mb is the absolute 
R-band SN la magnitude with X\ = 0, c = 0 and P = 0. 
The a, /3, <5 and Mb are the nuisance parameters which 
need to be fitted with the cosmological parameters. 

The JLA data set includes several low-redshift samples 
(z < 0.1), three season samples from the SDSS-II (0.05 < 
z < 0.4), and three-year data from SNLS (0.2 < z < 1). 
In total, it consists of 740 spectroscopically confirmed 
SNe la with high-quality light curves [18]. For the JLA 
data, the observed distance modulus is 

/ b LA = w”“ + a ■ xi - /3 ■ c - {M b + A M ), (2) 

and Am = 0 for M ga ,i < 10 10 M 0 , where M ga i is the 
mass of host galaxy. The Am denotes the correction of 
the absolute magnitude Mb for M ga i > 10 10 M 0 , and 
it is a nuisance parameter that would be fitted with the 
cosmological parameters. 

On the other hand, the theoretical distance modulus 
can be estimated by 

fi th {z) = 5 logi 0 UL ( 2 ) + 25. (3) 

Here Dl(z) = (1 + z) Dq{z) (in Mpc) is the luminos¬ 
ity distance from redshift z, where Dq is the comoving 
distance and is given by 

D c (z) = |U k r 1/2 sinn (^|H k | 1/2 ^ ’ ( 4 ) 

where f2 k is the cosmic curvature parameter, and 
sinn(x) = sinh(x), x, sin(a") for open, flat and closed cos¬ 
mic geometries, respectively. The H(z) = HqD,{z) 1 ^ 2 is 
the Hubble parameter where Hq = 100ft, kms _1 Mpc _1 , 
and f2(z) is expressed by 

Sl(z) =HM(l + 2 ) 3 + U DE (l + ^) 3 ( 1+ ^ + H k (l + z) 2 . (5) 

Here Om + ^de + H k = 1, and 12 m and Hde are the 
fractional energy densities of the matter and dark energy, 
respectively. The w is the static equation of state of dark 
energy. We have Ode = when w = — 1 for non-flat 
ACDM model, and f2 k = 0 for flat wCDM model. 

Therefore, the cosmological parameters are 12 m, Ua 
and ft in non-flat ACDM model, and 12 m, w and ft in 
flat wCDM model. Note that the Hubble constant Hq 
(or equivalently ft) has large degeneracy with the other 
cosmological parameters as shown in the H{z) defini¬ 
tion, and the SN la data actually cannot provide good 
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FIG. 1: The joint constraints on Dm vs. Da for non-flat ACDM model {left), and Hm vs. w for flat wCDM model {right). 
The constraint results from both of Union2.1 and JLA data sets are shown here. The solid filled contours and dashed contours 
are the results without and with systematic errors, respectively. The dots are the MCMC chain points used to illustrate the 
contour maps for the case without systematic errors. 


constraint on it. However, to avoid introducing priors 
and affecting the fitting results, we set Hq (or h) as a 
free parameter instead of fixing it to be 70 kms _1 Mpc _1 

We adopt the y 2 distribution to estimate the likelihood 
function C oc exp(—y 2 /2), and we have 


2 [^(Pl) -M th 0bPc)] 2 

X =2^- —I -> (6) 

2=1 

where N c j is the number of SN la data. The p c is the 
cosmological parameter sets, and p c = (17m, D a ,/i) and 
(Omi w, h) for non-flat ACDM and flat wCDM model, re¬ 
spectively. The p; denotes the nuisance parameter sets 
of the light curve, and we have p; = {a, /3,S, Mb) and 
{a, /3, Am, Mb) for Union2.1 and JLA, respectively. The 
a is the error for each SN la measurement, which is given 
by cr 2 = of c + cr 2 xt + <r 2 ys . The cri c is the error of light- 
curve parameters, and cr ex t includes the uncertainties of 
host galaxy peculiar velocity 3 and gravitational lensing 
effect 4 . Here we also include systematic errors cr sys in 
our estimation. Note that cr sys is obtained by different 
methods for Union2.1 and JLA data, which is derived 
by setting the reduced y 2 to be unity for each SN la 
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We assume 300 and 150 km s 1 for Union2.1 and JLA data, 
respectively BUI- 

We take <7i ens = 0.0932: and 0.0552: for Union2.1 and JLA data, 
respectively [161118j1 . 


subsamplc in Union2.1 and by REML method in JLA 
[2 12! 01 ■ The full covariance matrix between distance 
modulus are used for Union2.1 data with systematic er¬ 
rors, and we adopt the covariance matrix C in Eq. © 
instead of cr 2 in this case. 

We use the Markov Chain Monte Carlo (MCMC) 
to constrain the free parameters (nuisance parameters 
of the light curve and cosmological parameters). The 
Metropolis-Hastings algorithm is employed to determine 
the probability of accepting the new chain points pH . 0] , 
and the proposal density matrix is evaluated by a Gaus¬ 
sian sampler with adaptive step size 0. We assume 
uniform prior distribution for all free parameters, and 
their ranges in the MCMC fitting process are set to be as 
follow: D m e (0,1), D A e (0,2), w £ (-3,0), h £ (0.5,1), 
a £ (0,3), /3 £ (0,5), M b £ (-21,-17), S £ (-0.3,0.3) 
and Am £ (—0.3, 0.1). These ranges are chose n by our 
MCMC pre-runs and the relevant results from [T(J and 
171 . We find the fitting results are not very sensitive to 
the widths of these ranges. We perform sixteen paral¬ 
lel chains and get about 10 5 points for each chain after 
the convergence is reached [24| . After performing burn- 
in process and thinning the chains, we merge all chains 
together and obtain about 10,000 points to illustrate the 
probability distribution function of the free parameters 

In Fig. [□ we show the MCMC constraint results of 
Dm vs. D a and Dm vs. w for non-flat ACDM and flat 
wCDM models. We find Union2.1 data give similar con¬ 
straints to JLA data without systematic errors, and the 
best-fit values of the parameters are well consistent in l-cr 
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FIG. 2: The marginalized probability distribution function 
(PDF) of h for non-flat ACDM and flat wCDM models 
with systematic errors. The solid and dotted curves are for 
Union2.1 and JLA data sets, respectively. As can be seen, 
the constraints are loose for both ACDM and wCDM models, 
and the PDFs are relatively flat over the prior range. The 
results are similar for the cases without systematic errors. 


for these two data sets. When considering the system¬ 
atic errors, the deviations of the best-fit values of Dm, 
Da and w between Union2.1 and JLA become larger but 
are still within 1 -ct CL. We also notice that the contours 
with systematic errors are comparable to the case with¬ 
out systematic errors for JLA data. This could be due to 
the fact that there is no off-diagonal components for the 
covariance matrix of systematic errors between different 
SNe in JLA data, and just an cr sys (~ 0.1) is added to 
each SN subsample as the systematic error fl8| . 

We also investigate the constraints on h (or equiva¬ 
lently Hq) as shown in Fig. [2] Although there are peak 
features in the PDFs of h, we find it is not well con¬ 
strained by SN la data only as expected. The PDFs 
are basically flat over the range from 0.5 to 1 for both 
Union2.1 and JLA data. Here we emphasize that this 
flexibility of the constraints on h can provide unbiased 
constraints on the other cosmological parameters with¬ 
out introducing strong prior effects by fixing h or the 
Hubble constant Hq. Hence it could give reliable con¬ 
straint results of Dm, Da and w from SN la data only as 
well as the effective volume distance Dy. 

In Table (T) we list the MCMC fitting results for 
Union2.1 and JLA data with and without systematic er¬ 
rors, including both cosmological parameters and light- 
curve parameters. If we compare our JLA constraints in 
Table H| with the results in table 10 of 01, we can see 
that our constraints of [Dm, a, /?, Mb] on flat u>CDM 


model is consistent with the results in Ref. 0] within 
1 -ct CL. Then we compare our Union2.1 results with ta¬ 
ble 6 and 7 in 0, and we also find good agreements 
between these two results. Although the best-fit values 
of Dm and w are larger than theirs for flat u>CDM model, 
they are still consistent in 1 -ct. Note that our constraint 
errors are basically larger than that in 0 ] and 0 ] , and 
this is due to that we set h as a free parameter in our 
constraints and thus have one more free parameters than 
theris. 

In our constraint results, the cosmological parameters 
from Union2.1 and JLA data are consistent with each 
other in 1 -ct confidence level for the corresponding mod¬ 
els, but the light-curve parameters are a bit different since 
they are using different parameterization for the correc¬ 
tion between Ads and host galaxy mass. In addition, we 
also show the best-fits and 1 -ct errors for h in Table Q] 
but we should keep it in mind that its PDFs are rela¬ 
tively flat and it is not well constrained by SN la data. 
It is similar for Mb and A m in JLA data, that the PDFs 
are flat over the whole range, and just small and smooth 
bumps are shown in the PDF profile. 


B. BAO data 


The baryon acoustic oscillation is a good measurement 
for the cosmic distance and evolution. The BAO scale 
can be derived by fitting the scale dilation factor a, which 
depends on Dy/r d, according to a fiducial cosmological 
model. The Dy is the volume weighted effective distance 
at redshift z, which can be expressed as 


Dy{z) = 


(l + zfDl(z) 


cz 

H{z) 


1 1/3 


(7) 


Where Da{z) = Dc(z)/(1 + z) is the angular diameter 
distance, and we can also define a radial scale Dh(z) = 
c/H(z). The Da and Dr denote the distances can be 
measured by the BAO modes that perpendicular to and 
along the line of sight, respectively. 

We use several BAO galaxy clustering observations 
from six-degree-field galaxy survey (6dFGS) (9|, Sloan 
Digital Skv Survey Data Release 7 and 11(SDSS DR7 and 
DR11) 0 , 0 , WiggleZ dark energy survey 0 , and the 
Lya forest measurements from Baryon Oscillation Spec¬ 
troscopic Data Release 11 (BOSS DR11) 0, 0. The 
rd/Dy(z) values derived from these surveys are listed in 
Table HU 

The 6dFGS sample is obtained from more than 70,000 
half-sky samples out to z = 0.15 with effective redshift 
z e ff = 0.106. The BAO signal is detected at 105 h -1 Mpc. 
The SDSS DR7 Main Galaxy Samples (MGS) we use con¬ 
tain 63,163 galaxies and cover 6813 deg 2 at z < 0.2, and 
the survey gives Dy(z e e = 0.15) = (664 ± 25)(rd/rd,fid) 
Mpc 0 ]. The SDSS-III DR11 samples include nearly 
one million galaxies and cover approximately 8500 deg 2 
at 0.2 < z < 0.7, and they give Dy(z = 0.32) = 












5 


TABLE I: The best-fit values and 1 -a errors of the cosmological and light-curve parameters from the MCMC fitting for Union2.1 
and JLA data. 


Data 

Union2.1 

JLA 

Model 

non-flat ACDM 

flat wCDM 

non-flat ACDM 

flat wCDM 

Error 

no sys 

sys 

no sys 

sys 

no sys 

sys 

no sys 

sys 

9m 

Da 

w 

h 

a 

P 

Mg 

5 

Am 

0.26 ± 0.08 
0.65 ±0.13 

n 09 + 0.07 

U.M2_Q 38 
0.11 ±0.01 

2.31 ±0.05 
—18.70/q gg 
—0.03 ± 0.03 

0-28^ 

0-561q24 

0.77t° 0 f s 
0.10±0.01 
2.29 ±0.06 
-19.17±i;®® 

-0.14+S;£ 

0.28ZS:io 

-0.90t°;“ 
0-9lio,38 
0.11 ± 0.01 

9 

^.ou_o.o5 

— 18.70/q 
-0.03 ± 0.03 

n Qc;-l- U iU 

U.OO — 0.16 

—0 74 +°- 18 
u - '^-0.40 

0.75 ± 0.22 

0.10 ± 0.01 

9 9Q + 0 - 06 
^•^-0.05 

-19.62t;-« 
— 0.03 /q 24 

0.24 ± 0.08 
0.69 ±0.12 

0.88to'.37 
0.13 ±0.01 
3.13 ±0.08 
-18.51 

0.04l°;°« 

0.25 ±0.09 

0 . 66 ±g;ii 

0.57/“;" 
0.13 ±0.01 
2.96 ±0.09 
—18.50/“/® 

— 0.04/q24 

0.28/q4q 

—o. 92 / 5/22 

0 56 +0 ' 41 

U.OD_o 05 

0.13 ±0.01 
3.12 ±0.08 
—18.49/^24 

0.05/°;°| 

0 . 30/010 

— 0 .94/o 24 

0 72+ 0 - 27 
u - 1 z - 0.20 

0.13 ±0.01 
2.96 ± 0.09 

1 O ro+0.18 
IO. O 0 _o. 99 

—0 17+ 0 - 24 
U - 1 1 - 0.12 


TABLE II: The BAO data used in this work. 


Redshift 

r d /Dy(z) 

Data set 

0.1 

0.336±0.015 

6 dF [9] 

0.15 

0.2239±0.0084 a 

SDSS DR7 [26J 

0.32 

0.1181±0.0023° 

SDSS-III DR11 [10] 

0.57 

0.0726±0.0007 a 

SDSS-III DR11 [10] 

0.44 

0.0870±0.0042 

WiggleZ [27] 

0.60 

0.0672±0.0031 

WiggleZ [27] 

0.73 

0.0593±0.0020 

WiggleZ [27] 

2.34 

0.0320±0.0013 b 

SDSS-III DR11 [11] 

2.36 

0.0329±0.0009 b 

SDSS-III DR11 [28] 


“ The SDSS values here have been inverted from the 
published values of Dy(z)/rd, see the details in the text. 
b The BOSS values here are estimated from D\{z)/rd and 
Dn(z)/rd in the relevant references, and the details of 
estimation can be found in the text. 


are the square roots of the the diagonal values of C ma¬ 
trix. The two correlation coefficients are p = 0.54 and 
0.36 respectively, and the 1 st and 3 rd redshift bins are 
uncorrelated. 

Beside the galaxy surveys, we also include the measure¬ 
ments of the Lya forest results from DR11 of BOSS from 
SDSS-III at redshifts as high as z = 2.34 and 2.36 |ilTl.l28(. 
In Ref |l i|, the flux auto-correlation of Lya forest is used 
to detect the BAO features. There are 137, 526 quasars 
in redshift range 2.1 < 2 < 3.5 from DR11 of BOSS from 
SDSS-III are included to determine the position of BAO 
peak, and they find D\(z = 2.34)/r d = 11.28 ± 0.65 
and D h {z = 2.34)/r d = 9.18 ± 0.28. In Ref. H, the 
cross-correlation of quasars with the Lya forest absorp¬ 
tion is measured. The BAO scale is then given by Da(z = 
2.36)/r d = 10.8 ± 0.4, and D H (z = 2.36)/r d = 9.0 ± 0.3. 
Then, we use Eq. 0 to calculate r d /Dv and estimate the 
errors from the values and errors of Da/aj and -Dh/aj. 


(1264 ± 25Mpc)(r d /r d , fid ) and D y {z = 0.57) = (2056 ± 
20Mpc)(r d /r dj fi d ) 5 , where r dj fi d = 149.28 Mpc [l(|. To 
unify the measured BAO quantity as r^/Dy(z), we make 
conversion for the results of SDSS DR7 and DR11 and 
show them in Table HT1 The WiggleZ Dark Energy Survey 
samples contain about 200,000 redshifts of UV-selected 
galaxies, covering of order 1000 deg 2 of equatorial sky, 
and provides the BAO measurements at three redshifts 
M- In Jl^l, it is shown the 3x3 covariance matrix of 
Dv(rf d /r s ) for the three redshift points, and we convert 
this to the 3x3 covariance matrix of (r d /Dy) as 


C. CMB data 

Before the epoch of recombination, the baryons are 
coupled with photons by the large density of free elec¬ 
trons, and the photon-baryon fluid propagates as acous¬ 
tic waves. We can define the comoving sound horizon of 
this wave that can reach before the baryon decoupling or 
the end of baryon drag epoch. This sound horizon de¬ 
pends only on the physics of the early Universe, which is 
estimated by 


C = 


( 17.72 6.9271 0 \ 

6.9217 9.2720 2.2243 
\ 0 2.2243 4.1173 ) 


x 1CT 6 . 


( 8 ) 


Thus, the quoted errors of three WiggleZ data in Tabic ITT1 


5 This sample is the most precise BAO constraint ever obtained 

from galaxy survey. 


rd(zd) 



da . . 

a 2 H(a)(l + (3fl b /4D 7 )a) 1 / 2 ’ ( ’ 



where 2 d is the redshift of the drag epoch, c s is sound 
speed, t s denotes the epoch of last scattering, a is the 
scalar factor and a d = 1/(1+ 2 d ), and fib and are 
the fractional energy densities of baryon and radiation, 
respectively. 
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FIG. 3: The comparison of rd/Dv(z) data from BAO surveys 
with the prediction from Planck 2015 best-fitting cosmologi¬ 
cal parameters and l-cr errors (gray dashed line and region). 
The residual corresponding to the prediction of Planck 2015 
is also shown. 


The most state-of-the-art CMB observational data is 
the data from Planck 2015 results. Here we directly use 
the constraint on r d , i.e. the comoving sound horizon 
at the drag epoch, from the Planck 2015 results from 
cosmological parameters J2] 

r d /Mpc = 147.27 ± 0.31 {la CL). (10) 

This value is derived from the combined constraint of 
Planck TT, TE, EE power spectrum [||. We note that 
the r d is almost the same for ACDM and u>CDM models 
and the discrepancy is less than 0.15% [I, H|, since the 
dark energy is not important to affect the evolution of 
the Universe at early time. Therefore, we use the same 
value of r d given by Eq. (fllil) for both ACDM and w;CDM 
models in our following estimation. 


III. COMPARISON 

In this section, we compare and check the consistency 
between the results of the SN la, BAO and CMB mea¬ 
surements. We first check the consistency between the 
SN data and CMB measurements from Planck 2015 by 
comparing their x 2 and probability to exceed (PTE). 
Next, we compare the BAO results (r d /Dv) BAO with 
CMB results (r d /Dv) CMB - At last, we perform the com¬ 
parison using the quantity r™ B /Dy N with (r d /D v ) BAO 
and (r d /Z)v) CMB for all three kinds of observations. 

In Table Hill we show the results of x 2 and the PTE 
using Planck 2015 results (flat ACDM) as hypothetic 


model for Union2.1 and JLA data sets. The PTE is given 
by 


i r°° 

PTE = / f d/2-1 e -t/2 dt, (11) 

2 d / 2 T{d/2) J x 2 V ' 

where d is the number of degree of freedom, x 2 is from the 
model fitting and T(x) is the Gamma function. For com¬ 
parison, the minimum x 2 and corresponding PTE from 
our MCMC fitting for non-flat ACDM and flat u>CDM 
models are also shown. The degree of freedoms (dof) are 
572 and 732 for Union2.1 and JLA data, respectively 
(we have seven free parameters in all models). Note 
that since we have light-curve parameters fitted simul¬ 
taneously with cosmological parameters in the MCMC 
process, we use the results of Dm = 0.3156, Da = 1 — Dm 
and w = — 1 from Planck 2015 with the best-fit values 
of light-curve parameters (using Eq. (j6j>) to calculate the 
corresponding y 2 and PTE for Union2.1 and JLA data. 
We find that the x 2 and PTE of Planck 2015 are compa¬ 
rable with the minimum x 2 and PTE of non-flat ACDM 
and flat u>CDM models for both Union2.1 and JLA data 
with and without systematic errors. This indicates that 
the results of Planck 2015 and SN la data are in good 
agreements. Also, the PTEs for all models are large and 
greater than 0.48 which implies all these models are good 
to explain the SN la data sets. We note that some of the 
PTE values are quite high, especially for JLA data with 
systematic errors, which are close to 1. This implies that 
the models over-fit the data, and the errors of these data 
are probably overestimated. 

In Fig. [3j we show the comparison result of the BAO 
surveys and the predication of the best-fitting cosmo¬ 
logical parameters and l-cr errors from Planck 2015. 
The rd/Dy(z) values from the BAO surveys are listed 
in Table El The values of cosmological parameters from 
Planck 2015 we use to calculate rd/Dy{z) and errors 
are Dm = 0.3156 ± 0.0091, Da = 1 — Dm and H 0 = 
67.27 ± 0.66 kms _1 Mpc _1 |3|. As can be seen, the 
rd/Dy(z) of BAO measurements are consistent with the 
Planck 2015 best-fitting prediction in lcr CL, except for 
two data of SDSS-III DR11 at z = 0.32 and 0.57 and one 
data point of DR11 of BOSS from SDSS-III at z = 2.36, 
whose l-cr lower limits are a bit higher (within 2cr) than 
the Planck 2015 result. In order to estimate the com¬ 
bined significance of the consistency, we calculate the ef¬ 
fective x' 2 for the BAO data using the best-fitting cos¬ 
mological parameters from Planck 2015 as model pre¬ 
diction, which is given by 

N 

Xeff = [i r d/Dy)f MB - (r d /D v )F AO ] (Chao)*/ 

x [(r d /P v f B - (r d /Hv) BAO ] > (12) 

where N=9 is the number of BAO data, and Cbao is 
the full covariance matrix of nine BAO data. For 6dF, 
SDSS-DR7, SDSS-III DR11, their quoted errors are un¬ 
correlated, so they only have diagonal values in the Cbao 
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TABLE III: The x 2 and PTE for Union2.1, JLA data sets with the results of Planck 15 (P15, flat ACDM) as hypothetic model. 
For comparison, we also show the minimum \ 2 and the corresponding PTE for non-flat ACDM and flat wCDM models from 
MCMC fitting for Union2.1, JLA data sets. The degree of freedoms are 572 and 732 for Union2.1 and JLA data, respectively. 


Data 

Union2.1 

JLA 

Model 

non-flat ACDM 

flat tnCDM 

P15 (flat ACDM) 

non-flat ACDM 

flat uiCDM 

P15 (flat ACDM) 

Error 

no sys 

sys 

no sys 

sys 

no sys 

sys 

no sys 

sys 

no sys 

sys 

no sys 

sys 

X 2 

PTE (%) 

548.89 

74.94 

530.33 

89.31 

548.87 

74.96 

530.37 

89.29 

551.76 

72.10 

530.68 

89.11 

724.70 

56.90 

508.30 

99.99 

724.79 

56.81 

508.37 

99.99 

732.98 

48.28 

512.06 

99.99 


matrix ('Table HT1). But for WiggleZ survey, we use Eq. (0) 
as its covariance matrix. We find x 2 g = 7.93, the re¬ 
duced Xeff = 0.99 with dof=8 (we don’t have free pa¬ 
rameter here), and the PTE=0.44 with the Planck 2015 
results as the hypothetic model. This means the results 
of Planck 2015 can fit the BAO data very well. There¬ 
fore, generally speaking, the two kinds of observations 
are consistent with each other. 

Next, we perform our more robust test on cosmic evolu¬ 
tion by including SN la observations. For the numerator 
of r d /Dy (the scale of the sound horizon), since it is de¬ 
termined by the physics before 2 ~ 1100, we only use the 
CMB constraint from Planck TT+TE+EE power spec¬ 
trum which gives ru = 147.27 ± 0.31 Mpc |4j. For the 
denominator of r d /Dy, i.e. Dy(z) which is mostly de¬ 
termined by the low-redshift evolution of the Universe, 
we adopt the constraints given by supernova measure¬ 
ments from Union2.1 and JLA catalogs. Note that the 
Dy{z) from SN la catalogs is derived from our MCMC 
results. In order to avoid the statistical bias, we actually 
estimate 1 /Dy{z) for each MCMC chain point, which is 
composed of the cosmological parameters and light-curve 
parameters. We then derive the mean values and stan¬ 
dard deviations of t/Dy at different redshifts. Note that 
the mean values of 1 /Dy(z) are obtained by integrating 
over the full PDFs of l/Dy(z) from the MCMC results, 
which take account of the whole profiles of the PDFs. 

Then we combine the r d from Planck 2015 and the 
mean 1 /Dy(z) from SN la catalogs to construct the 
quantity r® MB /Dy N (z) 6 , and compare it to the results 
of [r d /Dy(z)] CKUi from Planck 2015 and BAO measure¬ 
ments [rd/Dv(^)] BA °. We show the comparison results 
in Fig. [I] For clarity, the residuals of r d /Dy{z) cor¬ 
responding to the predictions from Planck 2015 best¬ 
fitting results are shown. We derive rj MB /Dy N (^) for 
non-flat ACDM and flat ujCDM models for Union2.1 
and JLA data with and without systematic errors. We 
find the rj MB /Hy N (z) results are consistent with both 
[rd/Dv(z)] CMB and [rd/Dv(z)] BAO in Ict for all cases we 


Hereafter, for simplicity, we use the notation 1 //.l(A. 1/P| N and 
1 /D| n to denote the mean values of these quantities from the 
SN data, which are obtained by integrating over their full PDFs 
from the MCMC results. 


consider, and there is no strong evidence of the deriva¬ 
tion of standard ACDM cosmology evolution. The con¬ 
sistency of the non-flat ACDM model is similar to that 
of flat uiCDM model. Also, the values of rj MB /I?y N ( 2 ) 
are generally higher than [rd/-Dv(~)] BA °, which indicates 
the Dy(z) obtained from the SN la data is smaller than 
that given by the BAO measurements, especially at low 
redshifts (z < 0.3). 

In Table HVl we show the effective x 2 for rj MB /Dy N vs. 
{r d /Dy) CMB and (rd/Dy) BAO for different models and 
SN la data sets. The x 2 ff calculations here are similar to 
Eq. (fl^l) . when comparing r™ B /Dy N with (r d /Dy) BAO 

we use C = C , bao+C’sn,cmb: where Csn.cmb = Osn.cmb 
is the error of r d MB /.Dy N at each redshift. We find the 
effective x 2 is small and less than 4 for all cases, and the 
reduced x 2 g is less than 0.5. This indicates that all of the 
data sets are consistent with each other very well from 
the calculations of r d /Dy(z). Besides, using Eq. (fill) and 
Eq. (g2]) , we calculate the PTE for comparing r^ MB /Z)y N 
with (rd/.Dy) CMB with the results of Planck 2015 as the 
hypothetic model. Here, the x 2 i n Eq. m are derived 
from the xls by Eq. © with C = a 2 CUB + Csn.cmb 
where Oq MB is the CMB error from the result of Planck 
2015. We find the PTE are greater than 0.88 for all 
cases which means that (rd/-Dv) CMB is well consistent 
with Cj MB /J9y N . 

The BAO data at z = 2.34 and 2.36 from the auto¬ 
correlation of Lyct forest and cross-correlation of Lya and 
quasars are important for the comparison, since they are 
located at the median redshifts compared to the SN la 
surveys (z < 1.5) and CMB measurements. To make fur¬ 
ther comparison for these two data points, we directly 
use r d /D A (z) and r d /Dn(z) given by the Lya forest 
measurements 7 . We then derive the mean values and 
standard deviations of 1 /.Da(. 2) and 1/Dh(z) from our 
MCMC chain points for SN la data, and get the corre¬ 
sponding r^ MB /D| N ( 2 ) and r$ MB /D^(z) at 2 = 2.34 
and 2.36 with r d MB = 147.27±0.31 Mpc given by Planck 


' For uniformity, We inverse Da /v d and Dh A'a given by the Lya 
forest measurements to derive r d /D A (z = 2.34) = 0.0887 ± 
0.0051 and r d /D H (z = 2.34) = 0.1089 ± 0.0033, and r A /D A {z = 
2.36) = 0.0926 ±0.0034 and r d /D H (z = 2.36) = 0.1111 ± 0.0037. 























z z 


FIG. 4: Residuals of r&/Dv(z) by subtracting the predictions of best-fitting Planck 2015 results. The dark blue solid (dashed 
open) points with error bars denote rJ MB /Dv N (2) from Planck 2015 and SN la surveys without (with) systematic error. The 
BAO measurements are also shown with the same symbols as in Fig. [3] In order to make comparison with the BAO data 
clearly, we slightly shift the dark blue solid dots towards lower redshifts. Note that the results of rd MB /Z)y N (z) at z = 2.34 and 
2.36 are very similar with each other, and the two dark blue points with error bars are almost overlapped together at these two 
redshifts. Panel (a): Z?y N derived from Union2.1 data for non-flat ACDM model; Panel (b): P® derived from Union2.1 data 
for flat wCDM model; Panel (c): P® derived from JLA data for non-flat ACDM model; Panel (d): P>v N derived from JLA 
data for flat wCDM model. We find the rJ MB /Dv N («) from CMB and SN la observations is consistent with [rd/Dv( 2 )] BA ° 
from BAO measurements in lcr CL. 


2015. 

In Fig. [5j we show the comparison results of 
r C MB /^SN( z ) and r CMB /Z) SN ( ^ with [r d /D A (z)} BAO 

and [rd/.DH(£)] BAO from Lyct BAO measurements, re¬ 
spectively. Both of the cases with and without SN la sys¬ 


tematic errors are shown in dotted open and solid points 
with error bars. Similar to Fig. [H we compare the re¬ 
sults for both non-flat ACDM and wCDM models with 
Union2.1 and JLA data. We find the discrepancies are 
larger than lcr between Lya forest (purple inversed tri- 
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TABLE IV: The effective X' 2 f° r comparisons of the r d /Dy results for SN la, CMB and BAO data. The reduced X' 2 ff are also 
shown with dof=8. Besides, we calculate the PTE with Planck 2015 results as hypothetic model for comparing r d MB /D// <> 
with (r<i/Dy) CMB . 


Data 

Union2.1 

JLA 


r OMB 

/Dy" & (r d /Dy) 

(JMJL5 


/Dy" & (r d /D v ) H ™ 


/Dy" & (r d /Dv) 

(JMtJ 


/Dy" & (ra/Ov, 

BAU 

Model 

non-flat ACDM 

flat toCDM 

non-flat ACDM 

flat wCDM 

non-flat ACDM 

flat wCDM 

non-flat ACDM 

flat toCDM 

Error 

no sys 

sys 

no sys 

sys 

no sys 

sys 

no sys 

sys 

no sys 

sys 

no sys 

sys 

no sys 

sys 

no sys 

sys 

'2 

Xeff 

3.576 

3.942 

3.252 

3.736 

3.414 

3.722 

3.142 

3.508 

2.497 

1.633 

1.011 

2.121 

2.448 

1.636 

1.074 

2.085 

Xred 

0.447 

0.493 

0.406 

0.467 

0.427 

0.465 

0.393 

0.439 

0.312 

0.204 

0.126 

0.265 

0.306 

0.205 

0.134 

0.261 

PTE(%) 

89.32 

86.23 

91.75 

88.01 

- 

- 

- 

- 

96.18 

99.02 

99.81 

97.70 

- 

- 

- 

- 


angles) and Planck 2015 best-fits (gray lines), except for 
the data r&jD^iz = 2.34). The discrepancies shown here 
are clearer than the rd/Dy results shown in Fig. El since 
the Vd/Dy is derived by combining rd/D\ and rd/Dn. 
For the comparison including the SN la data, the results 
at z = 2.34 and 2.36 are quite similar with each other. 
The consistency of the ACDM model is similar to that 
of wCDM model, and they are consistent with all cases 
of the Lya forest measurements and CMB results in la 
CL. We don’t find significant deviation between SN la, 
BAO and CMB data. 


IV. SUMMARY AND CONCLUSION 

In this paper, we propose a robust method of testing 
the consistency of the cosmic evolution for a given cos¬ 
mological model. Our spirit is that if a model can fully 
describe the the cosmic evolution throughout the history 
of the Universe, it should be able to fit the observational 
data from the high redshifts to low redshifts. We find 
that the measurement of BAO is a sensitive indicator for 
this comparison, since it depends on both rd, which is 
determined from physics before recombination, and Dy 
which relies on the evolution of the low-redshift Universe. 
Then we can compare rd/Dy from BAO with the ra ob¬ 
tained from the CMB observations and the Dy derived 
from the low-redshift SNe la measurements, and check 
the consistency for these observations given a cosmolog¬ 
ical model. 

We use the SN la data from Union2.1 and JLA data 
sets, which are currently the two largest SN la data sets. 
In order to get the Dy and avoid bias on the constraint 
results, we free all parameters in the ACDM and wCDM 
models, especially for Hq, as well as the parameters of 
the light-curve. In the fitting process, we adopt MCMC 
technique and illustrate the probability distribution func¬ 
tion for each free parameter, and derive the mean values 
and standard deviations of l/Dy N at different redshifts. 
We use the BAO measurements from 6dfGS, SDSS DR7, 
SDSS DR11, WiggleZ and DR11 of SDSS from SDSS-III 
that the highest redshift can reach z = 2.36, and unify 
the BAO results to be [rd/-Dv] BA ° for all BAO data. The 
r CMB taken from Planck 2015 results, and we use the 


same rj MB for both ACDM and wCDM models, since it 
is quite similar for different dark energy models. 

We find that the observed BAO value is well consistent 
with the derived quantity of ra MB /Dy N over the redshifts 
from 0.1 to 2.36 for both ACDM and wCDM models. 
The consistency of non-flat ACDM model is similar to 
the flat wCDM model for all cases we consider. We also 
note that the ra MB /Dy N at z < 0.3 is higher than the 
[rd/Dy] BAO for both ACDM and wCDM models. This 
is probably because that the derived Dy N tends to be 
smaller than that from BAO at low redshifts. In order to 
further check the consistency of the BAO data at median 
redshifts (z = 2.34 and 2.36) provided by the Lya forest 
measurements, we also compare [rd/Dp/\ for trans¬ 
verse direction and [rd/DH] BA ° for radial direction with 
r d MB /-^A N an d r d MB /-^H N i respectively. We find that 
the r// MB /D^ for both of ACDM and icCDM models are 
in good agreements with [rd/DA] BA ° measurements in 1- 
cr CL. at z = 2.34 and 2.36. These results indicate that 
the SN la, BAO and CMB data are in good agreements, 
and there is no significant deviation from the standard 
ACDM model. Also, we should note that the uncertain¬ 
ties of the SN la data are probably overestimated in some 
cases as shown by the PTE values (e.g. see Table nnj. 
This implies that there can be potential deviations from 
the standard ACDM model hidden in the overestimated 
uncertainties, which need further confirmation by more 
accurate measurements from future SN surveys. 

In the future, the observations of 21-cm intensity map¬ 
ping will provide more BAO data across all different red¬ 
shifts. The BINGO telescope (BAO as Integrated Neu¬ 
tral Gas Observation) [H,[30] will cover 0.13 < z < 0.48, 
CHIME surve y ( Canadian Hydrogen Intensity Mapping 
Experiment) [3l| will cover 0.8 < z < 2.5, and future 
SKA (Square Kilometer Array) phase-1 will cover cosmic 
evolution up to redshift 2. The supernovae surveys, such 
as LSST (Large Synoptic Survey Telescope) which would 
provide a few hundred thousand or more SNe per year 
[33j, could greatly enhance the measurements of the evo¬ 
lution of the low-redshift Universe. All of these future 
observations will offer more accurate and reliable mea¬ 
surements on the cosmic evolution, and provide stricter 
constraints on the cosmological models. 
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FIG. 5: The comparison between the direct measurements of t^/Da and r^/Da from Lya forest and predicted values from 
SN la data sets and Planck 2015 at z = 2.34 and 2.36 with (dotted open circles) and without (solid points) SN la systematic 
errors. The grey vertical line is the prediction from Planck 2015 best-fitting cosmological parameters. We find the discrepancies 
are larger than la between Lya forest (purple inversed triangles) and Planck 2015 best-fits (gray lines), except for the data 
ra/DA^z = 2.34). The consistency of the ACDM model is similar to that of rcCDM model, and they are consistent with the 
Lya forest and CMB results in la CL. 
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